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ABSTRACT 


Solution algorithms are presented for the vertex m-center and the 
absolute m-center problem. Both algorithms use partitioning techniques. 
The algorithms use special properties of the max-min node to test for 
optimality. The vertex m-center algorithm establishes an order among 
all partitions of a graph according to the smallest vertex m-radius each 
partition can have. It then directs one to calculate the vertex m-radii 
only for those partitions which can provide a minimal vertex m-radius. 
The absolute m-center algorithm establishes an initial solution which 
may not be optimal. Other partitions are then tested against this 
solution to determine whether or not they provide a better solution. A 
point is reached at which no untested partition can improve the extant 


solution and the algorithm terminates. 
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le INTRO DUC HOM 


This thesis is concerned with the m-center problem of graph 
theory. It is desired to find m points on a graph which are closest in 
some sense to the nodes of the graph. If the m points are constrained 
to be nodes they are referred to as the vertex m-centers of the graph. 
If the m points are allowed to be anywhere on the graph, they are 
called absolute m-centers. Reed [3] presents a comprehensive summary 
of previous work done to formulate and solve the problem in various 
contexts. 

Analytical definitions for the vertex m-centers and absolute 
m-centers are given in section II, together with supporting definitions 
and notational conventions. 

Section III presents a solution algorithm for the vertex m-center 
problem which can be used if mis greater than one. The algorithm is 
Original in detail, to the best of the author's knowledge, although its 
general concept is not [1,2]. 

Section 1V contains the author's main contribution. A series of 
theorems is presented which leads to an absolute m-center algorithm. 
The algorithm guarantees that an optimal solution will be found and 
works directly from the vertex m-center algorithm delineated in section 
III. It is also shown that the absolute m-center algorithm of 


Rosenthal and Smith [4] does not provide an optimal answer in general. 








Section V summarizes the thesis and makes suggestions for 


further research. 








II. DEFINITIONS AND NOTATIONS 


Unless stated otherwise, the context of any definition is the 
connected graph G(N,A), where N is a non-empty set of nodes 
ivertices) and A is the set of undirected arcs incident to the nodes in 
feeelet X besthe set of all points on GiN,A). 

The two-tuple (i,j) denotes the shortest path in G(N,A) between 
points i and j, where the length of a path is the sum of the magnitudes 
piesociated with the arcs comprising the path. Let d(i,j) represent the 


length of (i,j). For any intermediate point k on the path between i and j, 
= 4 Ka) 

and 
Gh Slips) eel UH) mcmivel Vggagy) = 


DELINITION I: 








sf 
ij 


The distance matrix of G(N,A) is the square matrix 


where 


d.. =d(i,j). 
ij 


Let S be the set of all spanning trees of G(N,A). For Sy an 


element of S, let (i,j), be the shortest path between nodes i and j in 


k 


the spanning tree Si and let d(i,j), be its length. In general, 


kK 


Siyly SS hen) (1) 


ge 
This is true since removing arcs from G(N,A) can, at best, leave (i,j) 


unchanged. 








BELEINITION 2: 
a. The vertex center of a graph G(N,A) is that vertex 
@injme N, such that 


max d(i,C(v)) = min max d(i,j), i,j @€ N. 
i j i 


low The vertex radius of a graph G(N,A) is R(v), where 


R(v) = max d(i,C(v)), ie N. 
i 


DEEINITION 3: 
The central path of G(N,A) is the path TT) such that 


Voce NO ere. 


d(I,J),, = min max d({i, j) 7 


k i,j ss 

Rosenthal and Smith [4] were the originators of the central path 
concept. Their definition of the central path is subject to awkwardness 
caused by the possible existence of loops or circuits in a general graph. 
The author chose to define the central path in terms of spanning trees so 
as to avoid that difficulty. 

The spanning tree from which the central path is derived need not 
be unique and has several important properties. Let (i,j)* denote the 
path between node i and node j in the central-path spanning tree, and 


let d(i,j)* be its length. Because of (1), it follows that 


SUG ep) aan (e) Vr) ee (2) 








DEFINITION 4: 
a. The central-path vertex center of a graph G(N,A) is that 
vertex C(v,cp) € N, such that 


max d(i,C(v,cp))* = min max d(i,j)*, i,jé N. 
I j i 


be The central-path vertex radius of a graph G(N,A) is 
R(v,cp), where 


R(v,cp) = max d(i,C(vcp))*, ie N. 
i 


BHbOREM 1: 
For a graph G(N,A), 
R(v) <R(v,cp). 
Proof: It has been shown that 
ati }) Sd aes 
therefore 


max d(i,j) < max d(i,j)*, i,j N. 
i i 


It then follows that 


min max d(i,j) <min max d(i,j)*, i,j N. 
j i j i 


therefore 
R(v) <R(v,cp). 
EFINITION 5: 
a. The absolute center of the graph G(N,A) is that point 
C(ab) X such that 


max d(i,C(ab)) = min max (i,x), i@ N. 
i x €X i 
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b. The absolute radius of the graph G(N,A) is R(ab), where 


R(ab) = max d(i,C(ab)), i EN. 
iL 


THEOREM 2 [4]: 

The absolute center of the graph G(N,A) must lie at the midpoint 
of the central path of G(N,A). 

If C(ab,cp) denotes the central-path spanning-tree, absolute 
center of G(N,A) and R(ab,cp) represents the central-path spanning-tree, 
absolute radius, then Theorem 2 implies that 

C(ab) = C(ab,cp) 
and 
R(ab) = R(ab,cp). 
MEEINITION 6: 
The max-min node is that node V € N such that 


(V) = min d(i,V) = max min d(i,j) 


ifV j iFj 
It will be necessary at times to deal explicitly with distances 
from various nodes to the max-min node. Since a graph can have 
several max-min nodes the following convention will be used. Let 
V(A,B/1) denote the distance from node A to node B, where node Bis a 
max-min node of the graph and I is an ordering index. Then for any 
V(a,b/m) and V(c,d/n), if m <n, then 


V(a,b/m) <V(d,d/n) : 











DEFINITION 7: 

The graph G(N,A) is said to be partitioned if it is broken into a 
number of subgraphs. 

Let M be a set of m nodes chosen from G(N,A). The set M will 
be called a partition of G(N,A). For eachi EN, i ¢ M, place i in the 
set n associated with j ©M if 


d(i,j) = mind(i,k), k €M, 
k 


mas establishing 4 subgraph a ,A) of G(N,A) for node j. Each node 
j € M will be called a partition root of SINE 

The vertex center and vertex radius of GIN, ,A) will be C (vy) and 
R, (v) respectively where 


R.(v) = max d(i,C.(v))= min max dl(i,k). 
iéN, J k EN 1 EN, 


The absolute center and absolute radius of GIN ,,A) will be 
C (ab) and R, (ab) respectively where, for x, the set of points on 
GIN, ,A) ; 
R.(ab) = max d(i,C.(ab)) = min max d(i,x). 
i éNn. x€X. i én. 
) ) J 
et R, be the root length of GIN,,A) , defined by 


R. =maxd(i,j), iéN.. 
J a J 


Let Ruy be the partition radius of G(N,A) where 


Mme ie j}EM. 
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The absolute partition radius will be denoted R, ,(ab) where 


M 


Ry, (eb) = R, (ab), j€ M. (3) 


DEFINITION 8: 

The vertex m-center of the graph G(N,A) is that set of m nodes 
M* CN, such that Rus is minimized. That is, if R(m,v) denotes the 
vertex m-center radius, then 


R(m,v) = minR,.. 
MCN a 


DEFINITION 9: 


Let Xy be a set of m points on G(N,A) and define 


d(i,X,,.) =mind(i,x), i@N, x€X,.. 


M M 
x 


The absolute m-center of the graph G(N,A) is that set of m points 
* 
Xu CX, such that 


max d(i,X,*) =min max d(i,X,) : 
ién X= Kae & 


Let R(m,ab) denote the absolute m-center radius, then 


R(m,ab) = max d(i,X, *). 
ién 


Alternatively, the absolute m-center radius may be defined as 


R(m,ab) = min R,_ (ab). 
MCN 


Finally, let Sh be the spanning-tree from which the central-path 
of G(N,A) is chosen. Let C(ab) signify the absolute center of G(N,A). 


Given two nodes, i and j, j will be interior to i if and only if (i,C{ab)), 


i 








is identically (i,j ;C{ab))_. Conversely, j will be exterior to i if and 
only if iis interior to j. The node i will be an end-node of a path in 


S,. if and only if it is not interior to any other node inS 


K K" 


ue 








Wl. VERTEX M-CENTER ALGORITHM 


Reed [3] develops what might be called the brute force approach 
to solving for the vertex m-center of a graph. Using this method it is 
necessary to calculate Ruy for all possible partitions of G(N,A) and then 
to compare the Rag against each other to determine the minimum. 

The algorithm to be presented in this section appears to be more 
efficient. It makes use of the max-min node to pick those partitions 
of G(N,A) which might produce the minimum Ruy: For any partition of 
G(N,A) which does not have V as a partition root, the minimum distance 
from any root of the partition to V represents a lower bound on Rar’ 


Therefore the algorithm first finds R._ for those partitions in which V 


M 
is a root vertex. If it is then necessary to calculate Rus for any other 
partitions, the procedure tests partitions according to their increasing 


distance to V. Any partition having its distance to V greater than the 


R(v) is automatically rejected as a possible solution. 


AR THE ALGORITHM 
For a given graph G(N,A) with distance matrix inl ; 
1 Find those nodes which quality as V and note the value (v) : 
2. Determine aay for every partition having V as a partition 
root. If more than one node qualifies as V this step must be performed 


for each qualifying vertex. 


1s 








ce Find the smallest Ruy among those calculated in step D 
and test it in the following manner: 


a. IfminR,, < (V), go to step 4. 

lo Mbeseg ligt 12 = (V), go to step 5. 
—~M 

Celi mire \V)segenle step 7. 


M 


4. Set R(m,v) equal to min R Any partition among those 


ve 
determined in step 3 having Ruy = R(m,v) is a vertex m-center for the 
graph. Since any partition not identified in step 3 must have Ruy no 
smaller than (V) , the algorithm terminates. 


Br Set R(m,v) equal to min Rus: Any partition among those 


determined in step 3 having R 


vin R(m,v) is a vertex m-center for the 


graph. At this point a solution has been found and one or several 
partitions identified as vertex m-centers. Other partitions not having 
V as a root vertex may qualify as vertex m-centers also. 


6. Find those nodes j € N such that 


d(j,V) = (Vv). 
Again, this must be done for all those nodes qualifying as V. Determine 
Rag for all partitions having any such node j as a partition root. Any 


partition so determined whose R.__ is equal to the value of R(m,v) from 


M 
step 5 qualifies as a vertex m-center of the graph. Since all other 
partitions must have Rag greater than (V) , the algorithm terminates. 

ie It is necessary at this point to order the nodes of G(N,A) in 


terms of their distance to V. Suppose p and q were V nodes, then 


following the convention of definition 6, the list of values might be: 
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Vii, p/l), Vj,q/2), Vik, p/3), «+. 
The zero values should not be included in this ordering. It is possible 
for V(i, p/1) to equal V(j ,0/2). If this should happen, the remaining 
steps of the algorithm should be performed for all nodes which have 
an identical distance to V. 
On Start with the minimal V(A,B/1) . Calculate Rug for all 
those partitions having node A as a root vertex such that 


d(A,B) =min  d(wv,B). 
véeM 


If the minimum of the R,, so calculated is greater than the V(A,B/1) used 


M 
to begin this step, choose the next largest V(A,B/D) and calculate the 
appropriate Rar’ Continue doing this until 
; =V 
min Ruy (A, B/D) 


3) Set R(m,v) = min R,. as determined in step 8. Any partition 


M 
for which Ras was Calculated during any of the algorithm steps such that 
R . = R(m,v 
Ry (m,v) 
qualifies as a vertex m-center of the graph. Since any partition not 


considered already must have Ruy greater than the value V(A,B/1) of 


step 8, the solution terminates. 


B. FIRST EXAMPLE 

Consider the graph in Figure 1. Itis desired to find the vertex 
two-centers. Step | is to find V and (Vv) . Figure 2 shows the distance 
matrix for this graph. Each element in the row labeled "min," directly 


under the matrix, is the 


ts) 








Ines!) 


ij 
for the jth column. The underlined element in the min row is the max 
of the minimum values and identifies vertex 2 as V. The value (Vv) is 
seven. 

Figure 3 shows the partition distance matrix for those partitions 
having vertex 2 as a partition root. The column to the right of the 
matrix, labeled Ruy! has elements which are the maximum values from 
each corresponding row of the matrix. 

Since 

min R,, = ae 
the algorithm branches to step 5 resulting in R(2,v) = 7 and the 
partition fo j 5} being a vertex two-center for this graph. 

The algorithm cautions that other vertex two-centers may exist. 
This may or may not be important depending on the context of the 
problem. Looking at the distance matrix in Figure 2 it is seen that 
vertices 1 and 5 have distances to V equal to seven. Proceeding with 
step 6, Figure 4 shows the partition distance matrix for those partitions 
having nodes | and 5 as roots. From the Rus column it is seen that 


five other vertex two-centers exist for this graph and the algorithm 


terminates. 


e. SECOND EXAMPLE 
The graph for this example is shown in Figure 5, together with 


its associated distance matrix. The value corresponding to (V) is 


16 








underlined as before. Figure 6 shows the partition distance matrix for 


those partitions having node 5 as a root. It is seen that min R,. = (V) =5. 


M 
The other partitions which step 6 identifies are shown in Figure 7. Thus, 


R(2,v) is 5 and there are three vertex m-centers for this graph. 


D. DISCUSSION 

It is very difficult to make quantitative statements regarding the 
efficiency of this algorithm versus the "brute force" method. Too much 
depends on the characteristics of the graph under consideration. At 
worst, this algorithm might require one to calculate Rus for all possible 
partitions, just as the exhaustive approach does. At best, the algorithm 
can require as few as (N-1)! /(m-1)! (N-m)! calculations of Ray 


(where N is the number of nodes in N). This represents (m/N) of the 


(N)! /m! (N-m)! possible Ry calculations. 


e/ 








IV. THEvABSOLUTE ME@ENTRERER ORE 


The ideas presented in this section were motivated by the author's 
realization that neither the absolute m-center algorithm of Reed [3], 
nor the algorithm due to Rosenthal and Smith [4] guaranteed an optimal 
solution. Although Rosenthal and Smith state that their iterative 
technique will arrive at the optimal solution, this is not true in general 
as will be shown with a simple example. Reed did not claim that his 
algorithm determined an optimal answer in all cases and was unable to 
provide a test to determine when it failed. 

It is easily shown that the vertex radius of a graph iS an upper 
bound on the absolute radius of the same graph [3]. That is, 

R(v) >R(ab). 
To the author's knowledge, however, a lower bound other than zero has 
never been derived for R(ab). The theorems which follow provide such 
a lower bound and delineate its usefulness in attacking the absolute 
m-center problem. 
ies OREM 3: 

Let N be the number of nodes inN. If N 72 for G(N,A), then the 
C(v,cp) must be an interior node of the central path of G(N,A). If 
N <2, the C(v,cp) must be an end node of the central path. 

Proof: The case for N <2 is obvious. Suppose N >2 and the 
C(v,cp) is not an interior node of the central path. Let node q be the 
Civfepie fwo possibilities occur. 
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Case 1. Assume node q is an end node of the central path. Let 
k be any interior node of the central path and j be the other end node 
of the central path. Let m be the end node of any other path in the 
central-path, spanning tree. The following inequalities must hold: 
Sli ney) cael ult 
S10) eile =e ll < pen tay, 
GiGie) = = dime) += 
The first inequality establishes the fact that 


max dli,q)* = d(j,q)*. 
ieéNn 


The second and third inequalities show that any central path interior 
node will have a smaller maximum distance to nodes of the central 
path spanning tree than q does. This contradicts the assumption that 
node q is the C(v,cp) since if node q actually was the Civ,cp) then 


d(j,a)* = minumax dint. +1 € Ne 
f i 


Case 2. “As sume node q is either an end-node or an interior 
node on any path other than the central path. Let k be the first node 
interior to q which also is on the central path. Let m andn be the end 
nodes of the central path and j be the end node of any other path in the 
central-path, spanning tree. The following inequalities must hold: 

max (a(m,k)*, d(n,k)*] caiclie pine 
d(k,m)*< d(q,m)*, 
dike) < dia, a). 


This also contradicts the assumption that node q is the C(v,cp). 


Ng 








Since neither the assumptions of Case 1 nor those of Case 2 can 
be true, C(v,cp) must be an interior node of the central path. 
THEOREM 4: 

For a graph G(N,A), 

R(ab,cp) > R(v,cp). 

Proof: Let vertices q and k be end nodes of the central path of 
G(N,A). 

Case l. LiN<2, 

R(ab,cp) =3R(v,cp). 

Case Zandt NZ 

di(q, Wis =aciG et comm = diy icp), ki. 
Assume that 

MCS. WA CIS)) - Bae hr CIS) ahs) 
Now suppose that 

Rlab,cp) = zd(q,k)*< 2d(q,C(v,cp))* = sRiv,cp). 
Ii this is true, then 

Gi clic ena c Dike 
which is impossible. Therefore, 

R(ab,cp) >$R(v,cp) 
must be true. 
MaGOREM 5: 

For a graph G(N,A), let M bea partition, MCN. Let node j bea 
root of this partition. It must be true that 


R, (ab) > R,. 
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EnroO tee uo <2 the assertion is obviously true. Assume we ele Let 
nodes k and q be the end nodes of the central path of GIN A) . Define 
Ris such that 

Roe = max d(i,j)*, ee: 
Observe that 


RR .. 
al 


a: 


Case l. Assume j = Civ ,cp) Then 


R 


jx = max fel ee. Cll ill & Ri(v,cp). 
Therefore, from Theorem 4, 
R.(ab) >3R.,. 
J J 
Case 2. Assume j =k. Then 
R., = dG ,a)* = 2R, (ab), 


therefore 
R (ab) = ae : 
J J 
Case 3. Assume j is an interior node of the central path of 
G(N, »A) » but is not OIG) . Then 
R Sneee eli, Wer, Clinic! S cliente); 
therefore 
R.(ab) > 2R,, 
J J 
Case 4. Assume j is any node on a path other than the central 


path of GINA) ee ihen 


R 


j* 


=max ld(k,j)*, d(j,q)*)< d(k,q)*. 


Therefore 


Ze 








The proof is complete since cases 1 through 4 are exhaustive. 

This series of theorems is the foundation for the following 
absolute m-center algorithm. The algorithm first finds the minimum 
absolute partition radius among the vertex m-centers of the graph. The 
algorithm then uses the properties of the max-min node to find those 
partitions which could possibly improve the solution. If no such 


partitions exist the solution terminates. 


ies DHE ABSOLUIE M-CENITER ALG@ERITEM 

For the graph G(N,A) with distance matrix I, I : 

le Find all vertex m-centers of G(N,A) using the algorithm of 
section ITI. 


Ze. Determine R,,(ab) for each vertex m-center, applying 


M 
equation (3) and any absolute center algorithm which is convenient [3]. 
Ore Set R(m,ab) equal to the minimum value found in step 2. 
4, Determine the V(A,B/I) for G(N,A). For any V(A,B/I) such 
that 
V(A,B/1) < 2R(m,ab), 


calculate R,, for any partition M such that 


M 


d(A,B) = min d(v,B). 
v eM 


5. Calculate Bas (ab) for any partition from step 4 such that 
R <2R(m,ab). 
LS (m,ab) 
If at any time during these calculations a value of Ray (ab) is found such 


that 


Us 








Ray (ab) <“ R(m,ab) 


set R(m,ab) equal to this new lower value. Do not calculate R._ (ab) 


M 
for any M such that 
R.._ > 2R(m,ab). 
By (m,ab) 
6. When all the appropriate R, (ab) have been calculated the 


M 


Optimal solution is the extant R(m,ab). Any partition having R,, fab) 
equal to this value has a set of Ciel: which is an absolute m-center 
for this graph. Since any partition not already considered must have 


R,, (ab) greater than R(m,ab) the algorithm terminates. 


Be DISCUSSION AND EXAMPLE 

Consider the graph and distance matrix used for the second 
example in section III. This is a very simple graph, ... deceptively 
Simple. 

Rosenthal and Smith claim that their iterative technique must 
converge to an Optimal solution for the absolute m-centers of a graph. 
Their argument in support of this claim (on page 23 of [3}) hinges on the 
assumption an optimal solution has been found if it is not possible to 
improve the solution by moving some node from one subgraph to another. 
This is not true in all cases. 

Consider the graph in Figure 5. When the absolute m-center 
algorithm of [4] was applied to this graph to find the absolute two- 
centers, subgraphs as shown in Figure 8 resulted with R(2,ab) = 4. The 
algorithm presented above arrives at a better value for R(2,ab) for this 
graph. 


Z3 








It was found in the second example of section III that sets 
{2 : 5} : fi ,4} and {o 43 qualify as vertex two-centers for the graph of 
Figure 5. Forming the appropriate associations it is found that 41,5} 
produces the subgraphs depicted in Figure 8(a) while {p Sh and {2,4} 
produce the subgraphs of Figure 8(b), resulting in R(m,ab) = 4, as before. 
From the distance matrix it is seen that nodes l, 3, 4, and 5 have 
distances to V less than or equal to 2R(m,ab) = 8. Thus Rig should be 
found for {4,5} ,{3,5} , 2,5}, {1,5}, (4.33, &.2, &.3, G.2h, 
{1 3} and {3 , 2} : Rae and R,, (ab) have already been calculated for 
three of these partitions so that they need not be calculated again. 
Figure 9 shows the calculation of R,,. 


M 


Calculating R, (ab) for {4 5} , the subgraphs of Figure 10 are 


M 
produced, giving R, (ab) = 35. This is a better solution and immed- 
iately removes {1,5} ; {4 ,3} and {3,2} from consideration. Calcu- | 
lation of R,, (ab) for {3 5} ; {1,2} and {1,3} gives the same result as 
{4,5} . The algorithm terminates with R(m,ab) = 34. 

The solution process appears much more difficult and involved 
than it actually is. However, once the purpose of each step is under- 
stood, hand calculations can be made rapidly. 

The solution produced by the algorithm must be optimal. It 
starts with a solution minimally bounded from above by R(m,v) and then 
tests all m-tuples which could possibly improve the solution, using a 


lower bound criterion based on Theorems 4 and 5 to select those 


m-tuples. 
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V. CONCLUSION 


A. SUMMARY 

This thesis develops a vertex m-center algorithm and an absolute 
m-center algorithm form 71. The vertex m-center algorithm is a 
search technique using the properties of the max-min node to guide 
the search. This algorithm should substantially reduce the work 
necessary to find the vertex m-centers of a graph, especially when it 
is not necessary to find all such m-centers, because it avoids calculating 
the vertex radius of all possible m-tuples partitioning the graph. 

Several theorems are stated and proved which delineate a lower 
bound on a graph's absolute m-radius in terms of the vertex m-radius. 
This provides a test of optimality which, together with the demonstrated 
bounding properties of the max-min node, has been used to define a 
search technique for the absolute m-centers of a graph. 

This search technique must find the optimal answer since it looks 
at any partition which can possibly improve the initially derived answer. 
The algorithm avoids unnecessary calculations by using the lower bound 
developed for the Rai (ab) to eliminate partitions which could not possibly 


give a better solution. 


Be SUGGESTIONS FOR FUTURE RESEARCH 
The author believes that a solution process for the absolute 


center of a graph can be developed which uses the properties of the 


aS 








central-path end nodes to define an efficient algorithm. Such an 
algorithm has been developed by the author which successfully solves 
all examples given in [4]. It is assumed that the max-min node of a 
graph would always be an end-node of the central path, and proceeds 

to find the minimum, maximum length path through the graph having the 
max-min node as an end node. It has been shown that the process need 
not always work however, and more effort is necessary to establish the 
exact conditions under which it fails before attempting to use it 
generally. 

It seems inefficient to generate a distance matrix for a graph and 
to use this matrix for further m-center calculations. It should be 
possible to develop algorithms similar to minimal Spanning tree or 
shortest-route algorithms which work directly from the node-arc 
incidence relations and magnitudes to develop a partition having what- 
ever properties are desired. The work of Reed [3] leads in this direction 
and could be extended profitably. 

Finally, an effort should be made to extend the formulation 
context of the problem. Rosenthal and Smith [4] relate several possibilities 
and more should be available. One interesting idea would be to attribute 
probability densities to the arc magnitudes as is done in PERT and 


Similar techniques. 
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